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ABSTRACT 


The feasibility of using x-ray diffraction methods to measure re- 
sidual stresses in uranium and zirconium (Zircalloy-2) was investigated. 
A precision method was developed for the determination of diffraction 
peak positions and the precision associated therewith. The statistical 
tables of Fisher and Yates were used to determine what order polynomial 
provided the best least squares fit within the known precision of the 
observed data. It was found that a second order polynomial provided an 
adequate regression. With the aid of a desk calculator, less than five 
minutes'! calculation time is required to determine the peak position to 
a precision of + 0.01°. 


A General Electric XRD-3 Diffraction Unit was used in this investi- 
gation. The standard two-exposure method was employed in all stress 
measurements with the oblique position taken at 45°. Cobalt, copper, 
and chromium radiations were investigated in determining the best high= 
angle lines to be employed. The stress constants were determined using 
annealed specimens of uranium and Zircalloy-2 stressed in bending. The 
stress was measured with two strain gages placed one on either side of 
the irradiated area. 


The stress constant for uranium was determined Ко Бе 1308 + 110 
psi/0.01° shift in 420 for copper radiation on the(116)planes at 
29 = 158.3°. The stress constant for Zircalloy-2 was determined to be 
430 + 1 psi/0.01° shift inA2® for chromium radiation on the (10.4) 
planes at 20= 156.4”. 


The surface stress distribution adjacent to the butt weld of two 
flat plates of Zircalloy-2 was measured. It was found that the weld 
produced a significant effect upon the stress distribution adjacent to 
the weld. 


It is concluded that the expected precision as determined by 
statistical analysis can be attained using the method developed heréin. 
It is recommended that a round robin specimen be sent to several stress 
analysis laboratories to fully verify this fact and to establish confi- 
dence in this method of residual stress measurement. 


Thesis Supervisor: Robert E. Ogilvie 
Title: Assistant Professor of Metallurgy 
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I. INTRODUCTION 
Since the work of Lester and Aborn (1) in 1926, the use of x-ray 

diffraction techniques to determine residual stresses has been consider- 
ed feasible -- at least from the theoretical point of view.  Bragg's law 
as pertaining to diffraction of x-rays from a crystal may be stated as: 

пА = 2 d sin © (1) 
where n is an integer 

A is the wavelength of the radiation 

d is the spacing between atomic planes in the crystal 


O is the angle of incidence between the x-ray beam and 
the atomic plane. 


For plane stress as in the surface of a stressed specimen, classi- 


cal theory of elasticity predicts o 


o 
т d l+u sint Y (2) 
where 4d is the difference in the lattice spacing of the (hkl) 
planes parallel to the surface and at an angle V to 
the surface normal 
C is the stress 
E is the modulus of elasticity 
u is Poisson's ratio 


Combining the above two equations: 


=(E_ де = 
Be arene 4. oie (3) 


where K is the stress constant. 


1 For complete derivation see (2). 
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Equation (3) indicates that the stress in the surface of a speci- 
men is directly proportional to the difference between the peak position 
values of © as observed in the normal and oblique observations. The 
direction of this stress is defined by the intersection of the plane 
containing the x-ray beam in the two observations with the surface of 
the material. 

A great deal of work has been done in the past few years investi- 
gating mechanical failure of materials -- particularly in attempting to 
correlate the material strength and fatigue endurance to the physical 
and metallurgical characteristics of a material. Most investigators 
(3), (4),(5) believe the residual surface stresses to be an important 
parameter in this regard. 

Several methods have been devised to measure residual stresses. 

Of these, only the x-ray diffraction method is non-destructive. Another 
advantage is that only the stresses in the very surface fibres are 
measured as compared to the integrating effect of the various mechanical 
dissection techniques (6). However, the x-ray diffraction method has 
never become a generally accepted method of stress determination (6) be- 
Cause of the lack of precision in the repórted results. 

It is the firm belief of the authors that the method itself is 
quite reliable. , It is rather the improper experimental approach to the 
. problem by most investigators which has caused the difficulty. Since 


precision of the order of 0.01? is desired in measuring diffraction peak 





E 
positions, one must be extremely careful in the experimental procedure 
and, more important, in the treatment of the observed data. Precision 
is a mathematically defined concept (7) and one often overlooked or 
merely estimated roughly in research as well as in engineering work. 
With the exception of Ogilvie's method of least squares determination 
of the diffraction peak position (8), there is no evidence of an 
approach sufficiently careful to warrant the expected and assumed pre- 
Cision. It is little wonder that there has been a lack of precision 
and accuracy in the results reported using x-ray diffraction methods. 

When dealing with statistical data, one must apply the well es- 
tablished principles of statistics in order to obtain the maximum 
correct information from the observed data. The observed counting 
rate of diffracted x-rays is definitely statistical in nature as was 
verified by application of the Chi-squared test to a set of observed 
data. In experimental design it is appropriate to have the statistical 
precision limitations of the same order of magnitude or less than the 


inherent errors in the apparatus. 
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ІІ. PROCEDURE 

A precision method of diffraction peak position determination was 
developed. * The regression method by orthagonal polynomials (9), (10), 
(11) was applied to observed intensity data across the diffraction peaks 
chosen for subsequent stress evaluation. 65,536 counts were observed 
at each of seven equally spaced values of angle across the peak. Use of 
the Fisher and Yates tables (11) provides a quick method for determin- 
ing the order of Бо на. which should be fitted to best approximate 
the observed data," The criteria for fitting is the least squares fit 
which results in a mean square deviation of the same value as the 
variance of the observed data. It was found that a second order polynom- 
ial was the appropriate function to be used. ^ The direct method of peak 
determination as recommended by Ogilvie (8) was developed in terms of 
the orthagonal polynomials. This method in conjunction with the Fisher 
and Yates tables greatly simplifies this calculation for any number of 
points (up to 75 points;. Once the peak has been “counted,” less than 
five minutes’ calculation time is required by this method to determine 
the diffraction peak position from seven points. This least-squares 
method was used throughout this work. 

Uranium and Zircalloy-2 specimens were irradiated with x-radiation 


from cobalt, copper, and chromium in order to determine the appropriate 


2 Section C of Appendix. 
3 Section C of Appendix. 
а Section C of Appendix. 
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diffraction lines to be observed. As high a 2 6 value? as possible was 
chosen to get good precision. The reason for this may be seen as 
follows: 

From Equation (3), for constant stress and V 

д да tan 6 (4) 
which shows that a high-angle line is most desirable for best precision. 
The lines chosen weres 

Uranium Copper Radiation 20-7 158.3? (116) planes. 

Zircalloy-2 . Chromium Radiation 20-7 156.4° (10.4) planes. 

A General Electric XRD-3 diffraction unit was used throughout this 
investigation. The slit arrangement, filtering, high voltage and current 
values were varied to determine the optimum arrangement. In general it 
was desired to obtain as high a counting rate as possible to yield ade- 
quate precision in a reasonable time. No attempt was made to determine 
the background counting rate since no variation was to be expected over 
the region considered (0.30°). The final experimental arrangement 
chosen was: 

Filter -- none 


Divergence Slit == 3° with Solar baffles irradiating an 
area 6 mm x 1 cm on the specimens. 


Receiving Slit ==- 0.1? (positioned on the focus circle)® 


High Voltage == 47.5 kv 





5 20 the more commonly used parameter will be used hereafter. 
6 See (8) for complete discussion of the focusing conditions. 








Ga 
Current -- 16 ma 
Oblique Angle =- Y= 45? 
Collimating Receiving Slit == 0.3° 

It was found necessary to allow an hour's warm-up time for the 
machine with the x-ray tube on before a stable counting rate could be 
Aerea: ? This is most important for precise work. An appreciable 
shift in the apparent peak position will result if the average inten- 
sity is changing during the period of time required to "count" a peak. 
Also, the diffraction unit must be carefully aligned prior to precise 
stress evaluations. 

The stress calibration constants were determined by observing 
the change in A20 as the surface stress was varied known amounts. 
Specimens of uranium and Zircalloy-2 were vacuum annealed at 500? C 
for 8 hours to remove any residual stresses caused by fabrication. 
Туре A=7 SR4 strain gages were mounted on either side of the area to 
be irradiated. Each gage was read before and after each run, and the 
stress was calculated from the average strain value observed as 
measured with a type L Baldwin Strain Indicator ser J53470. A Sanborn 
Twin Visa continuous strain recording system was used to verify that 
the strain remained constant throughout a run. 

The specimens were stressed in a four-point bending device, 
thereby assuring a constant stress over the irradiated area. This was 


y Verified by Chi-squared test. 








verified experimentally by the strain gage measurements. The experi- 
mental setup is shown in Figure I, which shows the specimen in place 
and mounted on the goniometer of the diffraction unit. The alignment 
pedestal is shown in position for checking alignment of the specimen 
on the axis of the goniometer. Alignment of the specimen was checked 
before each run to within + 0.001 inches. 

The data for each peak was taken as follows. A slow scan of the 
peak was made. The apparent maximum point of the peak was determined 
to the nearest 0.05°. The time for 16,384 counts was recorded for each 
of seven points 0.05° apart centered about the estimated maximum value. 
Four separate counts across the peak were performed. This time sequence 
is more desirable than merely counting each point four times and then 
shifting to the next 2@value, as any transient effects tend to be 
eliminated. 

The geometrical correction factor due to any slight misalignment 
was taken to be the apparent A2 O0 value for a standard silicon specimen 
and a powdered tungsten specimen at approximately the same 20 angle as 
that used in the stress constant determination. Although this correction 
is relatively unimportant in evaluating the stress constant, K, it is 
of prime importance in absolute stress measurement. It was found 
appropriate to re-evaluate this constant for each experimental setup. 

In order to verify the precision method developed, each author evaluated 


the stress constant on a different specimen. 





ЕТСОКЕ Т 
Calibration Specimen Mounted 


on the XRD-3 Goniometer 











FIGURE I 


Calibration Specimen Mounted 


on the XRD-3 Goniometer 
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The weld-induced stress distribution adjacent to the tungsten 
electrode gas-shielded arc weld of two flat plates of Zircalloy-2 was 
measured. The 1? divergence slit was used in this stress analysis, 
limiting the area examined to 3 mm x 1 cm with the l cm side parallel 
to the weld. The welded specimen was mounted as shown in Figure II 
for stress evaluation. A stress pattern was determined as a function 


of distance from the center of the weld. 
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FIGURE II 


Welded Plate Mounted for Stress Analysis 
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FIGURE 11 


Welded Plate Mounted for Stress Analysis 
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III. RESULTS 

The results for uranium are plotted in Figure III. The stress 
calibration constant was evaluated as 1308 + 110 psi/0.01° shift in 
A 2е. 8 This value was determined by a weighted least-squares corre= 
lation with the predicted linear variation. The theoretical value 
was calculated to be 781 psi/0.01? shift of A2 a” The precision limits 
quoted for the stress constant as well as those indicated on the plotted 
points in the A20 direction were calculated from the known precision of 
the observed data. The precision limits for stress were estimated from 
the + 7.5 in/in observed precision in reading strain. The standard 
law of propagation of errors was employed in these calculations. 

The results for Zircalloy-2 are plotted in Figure IV. The stress 
calibration constant was evaluated as 430 + 1 psi/0.01° shift in A20 
using the same incon The theoretical value was calculated to be 
377 psi/0.01" shift A20 x 

Figure V shows the variation of stress in the welded plate. 

These stresses were measured in a direction normal to the direction of 


the weld along the centerline of the specimen. 





8 Е = 30x 10 psi. 
9 u = 0.28. 
10 Section B, Appendix. 


11 Е = 13.8x109 psi. 12 u = 0.32 
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Distance from weld centerline (cm) 


FIGURE V 


Plot of Welded Plate Surface Stress Distribution 











эте ыг 


IV. DISCUSSION OF RESULTS 

As shown in Figures III and IV, the results correlated quite well 
with the predicted linearity of stress as a function of 420. In each 
case the plotted curve does not pass exactly through the origin. The 
intercept of the curve with the ordinate axis indicates the residual 
stress in the specimen. This is to be expected, since it is extremely 
difficult to produce completely stress-free specimens. The deviation 
of the theoretical and experimentally determined values for the stress 
constants are probably due to the anisotropy of the material constants. 

The correlation between the observed and calculated precision is 
excellent. The data points, with but one exception, lie well within 
two calculated standard deviations of the least — MR straight 
lines. 

The precision limits of stress for uranium are greater than those 
for Zircalloy-2 because of the difference in the modulus of elasticity. 
The variation in the precision of A20 is due principally to peak broad- 
ening. This effect is included in the denominator term of the error 
equation. > 

It is interesting to note that even with the greater spread of 
data for the uranium-B specimen, the precision limits of the resultant 


stress constant were less than for the uranium-A specimen. Although 


13 Error equation is derived in Section C of the Appendix. 





= а 
this result is not to be expected in every case, the additional infor- 
mation provided by the extra point for the uranium-B specimen outweighed 
the comparative lack of precision for that specimen. 

The observed surface stress distribution in the welded plate 
(Figure V) indicates that the plate was not fully annealed prior to 
welding. This is substantiated by two factors. First, broadened 
diffraction peaks were observed in the regions away from the weld, the 
condition normally present in an unannealed plastically deformed material. 
Second, large compression stresses, which are characteristic of rolled 
plate, were observed far removed from the weld. The observed stress 
pattern is presumed to be due to the combined effects of fabrication and 


welding. 





з 1723 


V. CONCLUSIONS AND RECOMMENDATIONS 

It is concluded that with the proper experimental approach and 
Statistical analysis of the observed data, the x-ray diffraction method 
for measuring residual stresses is indeed a practicable one. This fact 
is especially emphasized by the excellent results herein obtained using 
diffraction peaks that would hardly be considered sharp peaks, as shown 
in Figure VI. The excellent correlation between the results from the 
A and B specimens confirms the validity of the method. 

It is recommended that this fact be firmly established by start- 
ing a round-robin specimen to various research and testing laboratories. 
This specimen should be accompanied by instructions for the method re- 
commended herein in determining the residual stresses in the specimen. 
An alternate procedure would be to re-evaluate the data and/or re-measure 
the stresses in previously circulated specimens (6). In this way the 
true merit of the x-ray diffraction method would soon become apparent. 

It is also recommended that further investigation be carried out 
in the fields ofs-- 


1. a comparative investigation of higher intensity lower 
angle lines. 


2. a Comparative investigation of photographic techniques. 


3. a complete stress analysis of the weld-induced stress 
in a properly annealed specimen. 


There are numerous possible applications of this technique in 
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c Ter 
the field of reactor design. Some of these would includes 


1. Evaluation of stresses produced by local hot spots in 
Zircalloy=2 cladding. 


2. Determination of the stresses induced by welding and 
fabrication techniques in fuel elements. 


3. Determination of the appropriate heat treatments to 
properly relieve residual stresses. 
The curve-fitting and error propagation techniques described 
herein are perfectly general. They can and should be applied in any 
experimental investigation, in order to properly establish the 


precision of the results. 
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REGRESSION AND ERROR ANALYSIS 
A. Regression Analysis 

Regression is the determination of the coefficients of an assumed 
mathematical expression purporting to represent the true means of an 
array of data. The selection of the regression equation may be left to 
the judgment of the analyst, in which case the simplest form producing 
adequate results is desired, or there may be physical laws which will 
govern the selection. Whatever form the regression takes it must be 
remembered that the methods presented here are only mathematical means 
to an end. 

Historically, the name, regression, originated with Galton, who 
first applied the method to determine how a population could remain in 
dynamic equilibrium if the sons inherited the characteristics of the 
father. He found that the sons! deviation from the mean stature was 
less than the fathers’ and consequently the sons "regressed" toward the 
mean. 

Consider an array of data with n columns and m rows such that one 
of the variables may be taken as independent. The first and probably 
most useful relationship between the variables would be that they were 
linearly related or nearly so. The questionsto be answered are: that if 
they are not exactly linearly related, what degree of approximation must 
be made to obtain such a relationship?, and, is the difference from 


linearity significant? Assuming X to be the independent variable, and 
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Y the dependent variable, the coefficients of the equation 
y! = ax + b (5) 
are to be determined. Yi; is the variate in the i th column and а 
j th row. The average value of the i th row is at 
> 


il 


зік 


У 2 = 


i Yij ° (6) 


=m 
= J 


J 


For reasons which will become obvious it is easier to determine the 


coefficients of 


y' =a (x - x) +b, (7) 
where x in the usual notation 
n 
p 
RESI 
zc ° (8) 


An appeal is now made to the Principle of Least Squares upon the assump- 
tion that it will lead to the most probable values of the coefficients. 


Therefore the minimum of 


n 2 
G (a,b) 2 (y Y) (9) 


is desired. Substitution of the assumed function yields 
n 2 
G (a,b) = 5 ÈA == b = aA (x, um 5) | о (10) 
1 


The minimizing conditions are 


n 
ASIS = à aA = 
22 lv. b = a (x, х)| 0 (11) 





n 
ЕЈ - 9 - 0, (12) 





п n 
>yY, =bn +a) (x, - X) (13) 
1 1 
1 E 
n n n 2 
2 (х, - X) у. DM (x, - x) +а (x - X) : (14) 
l 1 1 
Since by definition 
n 
У (х - ®) = 0 
1 1 
n 
a 
D = = 
xs Y (15) 
n 
20 7x) ^ 
and а + (16) 
n 2 | 
2. (x, - Xx) 
1 i 
with the resulting desired expression 
n 
y y = (m Leimen, , (17) 


n 
2 (х. - X) 


which is defined as the regression of y on x. 
With the constants of the equation now determined the question of 
the justification of a linear regression may be undertaken. The usual 


approach is to separate the deviation of any variate from the mean into 





= 26 = 
two parts, the deviation due to regression and the deviation from the 


regression. There follows 


n 2 n n 
= Y (у!- 7) +2) (y! - Y) (y, - y) +} (у = у!) 
19345 1 1 1 Д 
я 5 > 2 n E. n 2 
: Е) и: уа uU (y, - yn © 
since 
n 
LUE q, - Y) = 0. cu 


Thus the important result that the variance of the mean can be summed 
from the variance due to regression and the variance from the regression.* 
The variance from regression, (residual), is certainly a good measure of 
the justification of the linear regression as compared with some other 
assumed mathematical expression. Although there are no set rules it would 
seem inappropriate to select an expression of higher order if the average 


variance from regression is appreciably less than the average variance 





* It is entirely likely at this stage that an estimate of the statistical 
significance of the regression is desired. This may be obtained since 
there are now available two separate estimates of the same variance with 
n, and n, degrees of freedom. The "Z" or its equivalent "F" test or the 
better known Student's "t" test may be applied. 
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of the means of the variates. Mathematically this may be expressed as: 


2 


1 1 
5 (уы үл 72 (5 E 

2 1 = уа 
n 


п 


s |e 
ә ін 


(20) 


where 6 ы is the variance of the mean of the i th column in the array. 
If linear regression is found to be unsuitable, a logical exten- 
sion is to investigate polynomial regression. A straightforward approach 
would be to assume an equation of the type 
2 n 
y=b +bx+bx +... ÞXxX 
О 1 2 n 
and proceed as before by taking partial derivatives and applying the 


minimizing conditions. A better approach, as will be seen, is to define 


a set of orthagonal polynomials 


7 + + — 3 - г 
= а а х a x аж.+.а х 
упеге 
n 
Z 70 3 
2. 2. (22) 
n 
хао; (23) 
1 k 
n 
> 22 =0, n #k. (24) 
ı kn 
By this definition 
л = х 
1 (25) 
Б 2 n = 1 
27 = (хт ----) 
2 12 el 





RM 


3 Eno. 
я. Fe Es us) 
3 20 (27) 
2 2 
p. = 5 ЭПИ Зип) in?-9)| 
4 14 560 i (28) 


As an example consider fitting a second degree polynomial to the 
data. The assumed expression is 
у= А ҒА 2 + А 2. (29) 
O По 7 2 2 


The function to be minimized is: 





n 2 
G (Ap, Aj, A) - 2 (Уі = у а: (30) 
Upon substitution 
n 2 
С (А А А = = А = А 2 = А 2 
"mp T 2, EE оо 2’ | сш 
and differentiation 
А n 
2G =2 ОА ТАЙ -ла) (си | -о 2 
дА 2 - e l1 2 2 "M 
О 
n 
9G = 2 | = A -А2 -А2)(-2) F 0 
дА, 2 m О l1 2 2 1 = 
А n 
298, 22 ( = A -AZ -AZ2)(-2 | = 0. 
дА, 5 |(у, o 2 ща 
By equations (22), (23), and (24) 
n 
> (у -А) = 0; (35) 
1 1 
n 
A = l = y (36) 
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| 1 
n 
2. 
TRAG 
А, ^ n 2 (37) 
22 
ı 1 


Дт nee, (38) 


Further generalization leads to 


seen, 


2 У, 
P i 
n 


ма М5 


As mentioned previously the advantage of the method can be readily 


since the sums of the squares from regression ares 


Linear е А2 40 
217: та 

Очадга с: АУ 7 + А 77 (41) 
1 Te 2 24 

nth s AJZ + А52 + 2... ҒА 52 (42) 
MEI QM E n^ n'i 


The sum of the squares from regression can be easily determined and the 


effect of each additional stage of fitting evaluated. 


There still remains, however, the determination of the constants 
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of the orthagonal polynomials. The advantage of using orthagonal poly= 
nomials in the general case, when unequally spaced and/or unequally 
weighted values of the variates exist, is difficult to justify in so far 
i numerical labor is concerned; but for equally spaced and weighted 
values, the coefficients of the polynomials may be determined once and 
for all, and tabulated, for a cardinal number spacing. Tables of the 
polynomials for orders up to and including the fifth degree for up to 
and including seventy-five points have been developed by Fisher and 
Yates. These tables also include values for computing the sums of 


Squares necessary. 


B. Precision and Propagation of Errors 


Consider the problem of calculating the precision, or error, of a 
function made up of independently measured functions (i.e., empirical 
data) whose errors are known. Such a problem would be the calculation 
of the total resistance coefficient for a model, where the resistance, 
speed, and surface area cannot be exactly determined. 

The required assumptions ares 


а. Each quantity to be considered must be arrived 
at independently. 


b. The errors in the measured quantities must be 
small. 


с. Positive and negative errors must be equally 
probable. 


Let the true values of the measured quantities be X, S» ar etc. 
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These values are to be combined in some manner to give the desired 
function 


Х + Е (xi, Xj, X scent) (43) 


3? 
Let the true errors, which are unknown, of the measured quantities, be 
еі» е» еҙ ooo etc. The unknown error in X will be Е. Then 

X+E=f (x, + е,» X, t E x. + = E (44) 
By the required assumption (b) е,/х, is much less than one, anda 
Taylor's Series Expansion of equation (44) may be made with only the 
first order terms retained. 


Х + Е - Х- Е = (дх/дх.) e, + (дх/дх.) е, cece (45) 


By squaring both sides of equation (45) 


2 2 2 2 2 
E - (OX/ÓóÓx ) е + (9х/дж) е +... 
1 1 2 2 
ах/дх дх/д > 
+ (8X/ |) (дх/ х.) ее, , (46) 
where the terms contain the errors in the form е; апа SS Since (by 


assumption c) the errors are equally probable in a positive and negative 
sense, the cross product terms will tend to cancel out. Therefore, the 
desired result, which is commonly termed the Law of Propagation of Errors, 
15, 
Е дх/а ЕЕ 9х/9 =? 

= ( х/ х.) = + ( X/ ж) E Б 22222. (47) 
The error e of a value x is usually encountered in the form 

x + e 

This implies that there is a 50% probability that the true value of x 


lies between x plus e and x minus e. 
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С. Example 
Consider, as an example, the intensity data from an x-ray diffra- 


ction peak from a Zircalloy-2 specimen under stress. 


TABLE I. 


X-ray diffraction data 
Тарта Та 


36.56 35.42 34.87 34.27 34.83 36.22 36.49 
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In this array m= 4 
n= 7. 


Taking Z' values (Z' = ћ 2) from the Fisher and Yates tables for n = 7 
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TABLE II.' 


Coefficients of the orthagonal polynomials for n - 
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ТАВІЕ ІІІ. 


Calculation of variance 
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1335.9025 | -109.65 | 4182.75 +109.65 | = 36.55 
1264.5136 | - 71.12 | + 0 -248.92 | +142.24 
1225.7001 | - 35.01 | -105.03 + 35.01 | -175.05 
1197.8521 0 -138.44 +207.66 0 
1213.8256 -104.52 + 34.84 | +174.201 
1273.0624 + 0 ~249.76 | -142.72 
1329.3316 +182.30 +109.38 | + 36.46 
р Ганзен е еа оа 
у = 213. з5 53 
E 
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TABLE IV. 


Analysis of variance 


Degree of Fit Average 
0 ny” = 8836.6663 3.5216 ШШ 
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As discussed above, the residual, (variance from regression), is 
the required measure of goodness of fit. If a seventh degree equation 
had been chosen the residual would be zero, as the curve would pass 
exactly through each point. 

However the precision (i.e., the variance) to which each Уі 
value is known is easily calculable. For counting a constant number 


of counts (s) in time (t), the variance (52) is given bys 


52 -i (Е2 /в) (48) 


(The variance of the mean of four counts is one-fourth the variance of 
any one.) ** Therefore: 


for m= 4 


Е « 35 
s = 16384 
52 = „019 sec“ 


Thus it is seen that the second order equation is the appropriate 
regression to be employed. 

The information that is desired is not specifically the regression 
equation but its maximum value which may be determined as follows: 


уа = А + В! 21 + c'2, (in the notation of Fisher and Yates.) 
(49) 





** EVANS =-- pp. 756 and 798. 
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Determination of the precisions 


Assuming y у. +е 
1 qo 
ex 2 % 
x + Е = х + [ A | 
max ~ max “7 dy; i 
S 231. hy 
ee а is independent of 
> 2 zn. (is indepen y) 
li 
Assuming e. =e = 6 
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уз - T Суз Сти 2, ‚| 
^ 2 
(2 y,2),) 
2 TT = ng Z| бу а)? 22? + (туи Зан” 
- éy, (2 y,21) і2і ја ili 2i 
q = 
ы а е 
/ 
2 1:2 
Е = туђу у (2.1) + (2 у 2! j^ 2(2' ) | 
(Z yzi) 11 i li 2i 
i 2i 
= o 172 (2' j^ + х _ 2 (2! y * 
5 (жу, y? 11 max 2i (51) 


for n = 7 from Fisher and Yates 


h,-h -1 
1 2 
2 2 
à. (Z' )^ = 28 ; Liz!) = 84 
li 21 
т = 84 х 1 > 
28 x 1x 2 2 
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for the particular example 


x А20 84 l= „0176 
Max 17.06 2 2 


б = .137 


Lj 
| 


А 
TET E x 28 4.0003 x E 


00406 


In terms of spacing which is .05° 


х = „00088 + .000203° . 
тах = 


This calculated error due to counting statistics is negligible as compared 
to the geometrical precision of the diffraction unit. The instruction manual 
for the diffraction unit quotes a precision value of + 0.002°. Therefore, the 


final value of peak position is: 


156.050 4 9% 
max 


28 


20 156.651 + 0,0022 























